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O 1 ABSTRACT 

<n : 

Hipparcos parallax data for 40 contact binary stars of the W UMa-type 
(with eMy < 0.5) are used to derive a new, (B — V)-based absolute-magnitude 
calibration of the form My = My (log P, B — V ) . The calibration covers the ranges 
0.26 < (B — V) < 1.14, 0.24 < P < 1.15 day, and 1.4 < My < 6.1; it is based 
on a solution weighted by relative errors in the parallaxes (2.7% to 24%). Previous 
Q^ ■ calibrations have not been based on such a wide period and color space, and while they 

have been able to predict My with sufficient accuracy for systems closely following the 
well-known period-color relation, the new calibration should be able to give also good 
predictions for more exotic "outlying" contact binary systems. The main limitations 
of this calibration are the inadequate quality of the ground-based photometric data, 
and the restriction to the (B — V) index, which is more sensitive to metallicity effects 
than the (V — I) index; metallicities are, however, basically unknown for the local 
■ W UMa-type systems. 



1. INTRODUCTION 

Contact binary stars of the W UMa-type are - externally - simple structures. They are 
single objects with two mass centers and surfaces described by common equipotentials so that 
fewer parameters are needed to describe them than is the case for detached binaries: In place of 
two radii, one value of potential is needed, and in place of two surface temperatures, apparently 
one for the common surface suffices (Rucinski 1985, Rucinski 1993). The poorly understood 
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complexity of their internal structure is hidden from our view. The external simplicity and the 
observational property of the W UMa-type systems to appear close to the main sequence were the 
underlying reasons for establishing an absolute-magnitude calibration (Rucinski 1994 = CAL1) in 
terms of two observational quantities, the orbital period P and an intrinsic color, e.g. (B — V) or 
(V — I). Both quantities are correlated through the combined effects of similar geometry, Kepler's 
third law and main-sequence relationships, although the correlation is not perfect because of the 
evolutionary spread in sizes within the main sequence. 

The existence of a period - color - luminosity relation is somewhat reminiscent of that 
for pulsating stars. This should not be surprising if one realizes that the orbital periods for 
similar contact structures are in a simple relation to dynamical (or free-fall) time scales and 
that the radiating areas simply scale with the length of the period. A calibration of the form 
My = ap log P + ac C + const (where C is for any color index) can be also considered as an 
attempt to shift the stress from the color-derived temperature dependence of the normal main 
sequence (which presents operational difficulties, as it is affected by interstellar reddening) to 
a period dependence controlled by Kepler's third law. For the absolute magnitude calibration 
and derivation, the period can be considered essentially error-free, because it can be precisely 
determined by studying accumulated deviations of times of minimum light from a given ephemeris. 

As explained in CAL1, no mass-luminosity relation is involved in the calibration, just 
the simple geometrical relations. The most important omission consists of the neglect of the 
mass-ratio (q) term. The mass-ratio for newly discovered or faint systems is usually unknown as 
its determination requires spectroscopic observations at large telescopes, which are difficult to 
obtain for these typically faint and very short-period systems. Fortunately, the q-term possibly 
correlates with, or can be absorbed by, the two other terms. The other limitation of any luminosity 
calibration for contact binaries is that it is likely affected by the presence of spots on such very 
active stars. Therefore, the calibration applies to a somewhat poorly defined mean maximum 
brightness, where the averaging is over a spot-induced variability of poorly known temporal 
characteristics. 

The first calibration (CAL1) was based on 18 systems, mostly members of high galactic- 
latitude open clusters (to avoid Milky Way interlopers in more typical open clusters), and of 
visual binaries, as well as on 3 nearby systems with known trigonometric parallaxes, 44i Boo B, 
VW Cep and e CrA. Even for these three systems, the data were somewhat uncertain: The color of 
44i Boo B had not been measured with any degree of accuracy due to the light of the close visual 
companion, and large discrepancies existed between published parallax values for e CrA. From the 
beginning, there was a need for reliable trigonometric parallaxes for a large number of systems. 

CAL1 concentrated on the My = My (log P, B — V) relation. Further work related 
to an approximate evaluation of the metal abundance effects, based on systems in very 
metal-poor clusters (Rucinski 1995 = CAL2), and to a calibration of the relation in the form 
Mi = Mi (log P, V — I), to study the galactic distribution of systems discovered in the OGLE 
microlensing search (Rucinski 1997a = CAL3). In summary, the established relations were: 



My = -2.38 log P + 4.26 (B - V) + 0.28 - 0.3 [Fe/H] 
My = —4.43 log P + 3.63 (V - I)o — 0.31 — 0.12 [Fe/H] 



(1) 
(2) 
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Mj = —4.6 logP + 2.3 (V — I)o — 0.2 — 0.12 [Fe/H] (3) 

where the first two equations combine the results of CAL1 and CAL2, while the third combines 
the results of CAL2 and CAL3. The colors in Eqs. || and || are de-reddened and the (V — I) 
color is measured in the Kron - Cousins system. The (V — I)-based calibrations (Eqs. ^ and 
^) must be considered preliminary because of the lack of directly observed colors. Due to 
strong inter-parametric correlations, the coefficients in the relations have very large and highly 
non-Gaussian errors (see CAL1 and below, Section |3|). However, the absolute magnitudes can be 
predicted quite accurately. As was illustrated through Monte Carlo simulations, the preliminary 
Mj calibration (Eq. |||) could predict the absolute magnitudes to one-sigma (68.3% confidence) 
level of ±0.25 magnitude and to two-sigma (95.4% confidence) level of ±0.5 magnitude (see 
Figure 3 in Rucinski (1996)). At first, the period range was only 0.27 < P < 0.59 day, but was 
slightly extended to 0.63 day in CAL3. 

The extant versions of the luminosity calibration have been used by several authors, mostly for 
consistency checks on the membership of individual newly discovered systems in various clusters 
(Edmonds et al. 1996, Kaluzny et al. 1996a, Kaluzny et al. 1996b, Mazur et al. 1995, Rubinstein 
and Bailyn 1996, and Yan and Mateo 1994). A somewhat more ambitious application was its use 
for distance determinations within the "pencil-beam" search volume of the micro-lensing project 
OGLE (CAL3). 

Although the predictive power for individual cases is relatively low (conservatively estimated 
to be more accurate than 0.5 magnitude, taking into account the possibility of systematic errors), 
the various forms of the the calibration have a definite potential. Of particular importance is the 
high frequency of occurrence of the contact binaries: Disregarding the difference in population 
characteristics, the contact binaries are some 24,000 times more common than RR Lyr stars 
(CAL3) in the solar neighborhood. Thus, one can imagine their application to detailed galactic 
structure studies, in place or in addition to those based on RR Lyr stars. In addition, an extension 
of the calibration towards early-type contact systems may enhance their usefulness for more distant 
stellar systems. The calibrations may be also useful in obtaining a better estimate of the spatial 
density of contact systems. The apparent density derived from the OGLE sample (CAL3) of about 
one such system per 250 - 300 main sequence stars of similar spectral types (F-K) is about 3 to 
4 times higher than the density estimated from a sample of systems in the solar neighborhood 
(Duerbeck 1984). The discrepancy may be due to incompleteness of the contact binary sample 
below the brightness level of about V = 8. The complete sample of the F-K stars in the Hipparcos 
Input Catalogue (Turon et al. 1994) contains 128, 497, 1620 and 4561 stars of luminosity classes 
IV and V within the range of occurrence of the W UMa-type systems, 0.4 < (B — V) < 1.2 to 
magnitude limits V = 5, 6, 7, 8. To the same brightness limits, there are 3, 3, 6, and 14 contact 
binaries giving the inverse relative frequencies 43, 166, 270, 326. Thus, the frequency derived from 
the OGLE data is approximately confirmed, but with a large statistical uncertainty and with a 
trend indicating that discovery selection effects become progressively stronger for fainter systems. 
Once a larger sample of such faint (i.e. not accessible to direct parallax measurements) contact 
systems is identified, the calibration can be used to define a volume-limited subsample for an 
unbiased local density estimate. 
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Parallax data from the Hipparcos satellite (ESA 1997) have already been used to derive a 
new absolute-magnitude calibration for the W UMa-type systems (Rucinski and Duerbeck 1997 
= CAL4). The main motivation was the obvious need of having a solid check on the previous 
calibrations which had been based on inhomogeneous data from various sources. A sample defined 
by the errors in My smaller than 0.25 magnitude (eMy = 2.17^ < 0.25) consisted of 20 systems, 
but the system with the largest and most precisely measured parallax, 44i Boo B, still does not 
have a reliably measured color and had to be omitted. Only 2 systems were common to CAL1 and 
CAL4. The calibration directly comparable to CAL1 is CAL4a: 

My = -4.30 log P + 3.28 (B-V) + 0.04 (4) 

CAL4a turned out to show a small standard error (a = 0.17) and, despite its somewhat differing 
coefficients from those of CAL1, to give encouragingly similar predictions on My as compared 
to the previous calibrations. Monte Carlo simulations showed that, as far as random errors are 
considered, absolute magnitudes of contact binaries can be predicted with an accuracy of about 
0.1 magnitude, which is surprisingly good in view of the omission of the mass-ratio term and of 
the possibility of star spots. However, deviations of individual systems from the "fundamental 
plane" of the calibration indicated the possibility of systematic rather than random errors, and the 
problem of the applicability of the calibration over wider ranges of periods and colors remained. 
CAL4 was established within the period range 0.27 < P < 0.65 day, with most calibrating 
systems clustering around the main period-color relationship. The experimental addition of the 
long-period, red system V371 Cep (V5 in NGC 188) led to very different calibration coefficients. 
Since the two quantities of the calibration, period and color, are strongly (but not perfectly) 
correlated, the calibration is extremely sensitive to the inclusion of severely spotted or non-contact 
systems mimicking genuine W UMa-type systems. Realizing these limitations of CAL4, we 
consider in the present paper a larger sample, comprising also systems with poorer parallax data, 
with a hope that such a sample would cover a wider range of periods and colors, and thus lead to 
better averaging of individual peculiarities than in the case of the 19 systems used to establish 
CAL4. 



2. Hipparcos and auxiliary data 

Our present sample consists of 40 systems with parallax errors implying uncertainties 
eMy < 0.5 magnitude. Most systems are well-known, bright W UMa-type binaries. While details 
on the adopted photometric data for the sample systems are given in the next section and in the 
extensive Appendix, several systems included in or excluded from the sample require a special 
mention: 

44i Boo The color data for the contact binary which is the fainter component (B) of a close visual 
system are too uncertain to warrant its inclusion in the sample. 44i Boo B is the nearest 
contact binary with a parallax of 78.4 ±1.0 mas (milli-arcsec). 

BW Dra We included this system although the relatively uncertain parallax data (15.3 ± 5.2 mas) 
would formally exclude it from consideration. This system forms a visual binary with 
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another contact system, BV Dra. The better-quality parallax of BV Dra was assigned to 
BW Dra, too. 

AP Dor The poorly known contact binary HD 33474 was included in the sample although only 
a fragmentary light curve and a coarse estimate of the period existed until now (Eggen 
1980). The value of the period has been taken form the Hipparcos Catalogue (ESA 1997) 
where it is classified as a W UMa or possibly a RR Lyrae star. However, Eggen (1980) 
argued convincingly that AP Dor is a contact binary system. Its absolute magnitude also 
contradicts the RR Lyrae classification. 

BV Eri This binary was initially included in the list of contact binaries to be studied with Hipparcos 
(HIC 18080), but it has been excluded from our present study since the photometric analysis 
of Baade et al. (1983) has shown that it is a close detached binary. 

UZ Oct This contact system has period slightly longer than one day and formally does not belong 
to the class of W UMa-type systems. We included it in our calibration. It contributes little 
to the final determination, but its deviation from the calibration established for short-period 
systems may give us indications of the expected trends for long-period systems. 

ER Ori According to the photometric and spectroscopic analysis, ER Ori should have a trigonometric 
parallax between 5 and 7 mas, which would probably have lead to the inclusion into our 
sample. The Hipparcos parallax is, however, —6.8 ± 1.4 mas. Since ER Ori is a triple system 
with a third component at a distance of 0".l, found by speckle observations (Goecking et al. 
1994), it is likely that the orbital motion around the common center of mass has resulted in 
an unrealistic parallax value. 

In addition to the parallaxes and the periods, the input data for the calibration consist of 
the observed maximum magnitudes in V, the (B — V) colors at maximum brightness, and the 
reddening values E(B — V). No measured values of (V — I) were available for most systems, 
thus no separate My = My (log P, V — I) calibration was attempted (in most cases, the (V — I) 
data listed in the final Hipparcos Catalogue have not been measured, but were derived from 
{B — V)). This color index has become a standard one in CCD-based studies of stellar clusters 
and micro-lensing searches, and it is unfortunate that we presently do not have a good database 
to carry out a calibration. 

The V magnitude and (B — V) color index at maximum light were usually taken from the 
literature, but sometimes had to be inferred from Hipparcos magnitudes or from Stromgren 
photometry. Details on all systems are given in the Appendix to permit a quality assessment 
in each case. The adopted values of the data used to establish the new calibration are listed in 
Table |]. It should be stressed that some of these input data have gone through a much more 
careful scrutiny and thus supersede those adopted in CAL4. Although our main goal was to 
establish reliable values of the magnitudes and colors at light maxima, both Table |l] and the 
Appendix update and extend a very useful compilation of parameters of contact systems by 
Maceroni and Van't Veer (1996), which focussed on contact systems with modern photometric 
solutions. Our Table includes 14 systems not listed by them - in most cases because no detailed 
analysis has been carried out until now. 
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As indicated, the available V max data for some systems were poorly known or considered not 
trustworthy. In such cases, we gave a higher weight to independent estimates of 7 max , based on 
the Hipparcos data for the upper 5% of the brightness (the 95-percentile levels), Hp mSLX . Since the 
Hp photometric system slightly differs from V, a transformation to V magnitudes was necessary. 
In CAL4, we used a Hp to V transformation involving a small (V — I)-term which was estimated 
from the adopted (B — V) via the color - color calibrations of Bessell (1979). This two-step 
process has been simplified here by considering differences Hp — V directly as a function of 
(B — V). Reliable values of Hp — V at maximum light for 30 of our systems were used to establish 
two- and one-parameter relationships, Hp — V = 0.17(5 — V) + 0.03 and Hp — V = 0.22(5 — V). 
We adopted the latter, noting that the intercept should be zero by definition, but that in the 
present case it may be non-zero because of the averaging at the 95-percentile level in Hp. 

If no (B — V) index was available, it could sometimes be calculated from a measured (b — y) 
index, using relations established by Bessell (1979). Some (B — V) indices were taken from an 
unpublished photometric survey of southern contact binaries (Duerbeck 1997). In all instances, the 
data were compared with the mean data from the Tycho experiment tabulated in the Hipparcos 
Catalogue. 

The reddening E(B — V), although a relatively small quantity for the nearby systems 
considered here, can affect the calibration in a systematic way. It tended to be over-estimated in 
previous studies (Rucinski and Kaluzny 1981, Rucinski 1983) through the use of general distance 
and galactic-latitude scaling laws, which may be invalid for small distances. In fact, the mean 
value of the extinction of Ay = 1.8 mag/kpc (Whittet 1992) applies only to situations of good 
averaging of many individual interstellar clouds at large distances. For moderate distances below 
100 pc, and well within the galactic disk, the interstellar medium shows strong inhomogeneities, 
with some directions almost devoid of any absorption. The novel approach, adopted in CAL4 and 
in the present paper, uses values of neutral hydrogen column density, _Vhi, to estimate E{B — V). 
A large database of determinations for JVhi is available as an Internet hydrogen column density 
search tool from the Center for Extreme Ultraviolet Astrophysics in Berkeley, California^]; it 
is based on the data in Fruscione et al. (1994). On the basis of the studies by Bohlin et al. 
(1978) and Tinbergen (1982), the adopted relation between E(B — V) and Nui (cm -2 ) was: 
E{B - V) ~ 1.7 x 10~ 22 iVHi. A newer study, based on X-ray data (Predehl and Schmitt 1995), 
gives an almost identical relation with the coefficient of 1.8 x 10~ 22 , a difference which can be 
entirely neglected for the small values of E(B — V) encountered in here. The values of E{B — V) 
in Table [l] derived in the above way carry uncertainties of about 0.01 - 0.02 magnitude. On the 
average, they are much smaller than previously assumed for the field W UMa systems, but remain 
uncertain due to interpolation of the values of iVm over an irregular grid of stars, with relatively 
large spacings in angular separation and distance. 

The absolute magnitudes have been derived from My = V + 5 log7r — 3.1 E{B — V) — 10, 
where the parallax tt is in milli-arcsec units (0".001), and the observed magnitude at maximum 
light is in the Johnson U-band. 



4 http://www.cea. berkeley.edu/~science/html/sci_archive_tools_tools. html. 
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3. The new absolute-magnitude calibration 

The color - absolute-magnitude diagram for the systems used in this paper is shown in 
Figure [I]. Note the vertical spread of the points, which is mostly due to the evolution of some 
systems away from the main sequence and to un-accounted mass-ratio differences. Without this 
spread, the calibration could be simplified to just one independent variable. The period - color 
relation is shown in Figure ||. In these and the following figures, filled circles mark the systems 
with better parallax data (eMy < 0.25) whereas open circles mark the systems with larger relative 
errors in the parallaxes (0.25 < eMy < 0.5). There are almost exactly the same numbers of 
systems in each of these categories. 

The combined 3-dimensional period - color - luminosity relation is shown in Figure ||. The 
planar fit corresponds to the calibration 

My = a P{BV) log P + a BV (B - V) + a (BV) > (5) 

whose coefficients are listed in Table |2| The deviations from the fit are shown in Figure ||. As 
in the previous cases, the coefficients in Eq. |B| have errors which are very large and very strongly 
non-Gaussian. In CAL1 and in the following studies, the errors were estimated with a technique 
of "bootstrap re-sampling" , which utilizes a large number of solutions based on the observational 
data re-arranged randomly into several data sets of the same size by a random (and repetitious) 
selection of points. This technique permits to evaluate realistic errors for cases of their entirely 
unknown distributions. The application to our data is illustrated in Figure |B] which gives the 
spread in the coefficients obtained from 10,000 bootstrap re-sampling solutions. The one-sigma 
ranges in the diagrams can be identified with the extent of contours containing 68.3% of all 
cases. The essential statistics for the bootstrap re-sampling results is given in Table [| Note that 
the median and average values of the coefficients are closer than in the previous calibrations, 
indicating a better overall stability of the current solution. Since this calibration is based on a 
larger sample than all previous ones, Figure || no longer shows multiple "islands" of the solutions 
which indicated that the previous versions of the calibration were sensitive to the inclusion or 
rejection of individual systems. 

The relatively large uncertainties in the coefficients do not signify that the predicted values of 
absolute magnitudes are equally uncertain, as inter-parametric correlations result in cancelation of 
the contributing uncertainties. To determine uncertainties in the predicted absolute magnitudes, a 
Monte Carlo experiment was carried out. This experiment is based on the availability of the 10,000 
bootstrap solutions which were used for evaluation of the coefficient errors, as described above. 
We added one more level of randomization here by considering 100 random values of (B — V) 
for each value of the period between 0.2 and 1.0 days in steps of 0.01 day. This way, for each 
step in the period, 10 6 random combinations of the colors and coefficients were used. Whereas 
the distribution of the coefficients had been given by the bootstrap re-sampling experiments, the 
distribution of the colors was assumed flat within period-dependent ranges. These ranges increased 
linearly in width from A(B — V) = 0.3 at the period of 0.2 day to A(B — V) = 0.6 at the period of 
one day. The widening was meant to represent the observed trend in the period - color relation, 
with a much larger spread in colors for long-period system. This spread is usually interpreted as 
due to evolutionary effects which are more pronounced for long-period systems. The blue edges 
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of the respective color ranges were placed at the observed blue, short-period envelope (BSPE) 
of the period - color distribution. The existence and importance of the BSPE was discussed in 
CAL3 where an approximate fitting formula was determined: (V — I)bsp = 0.053 P . Here, 
we used the same BSP envelope, but transformed to the {B — V) color using the Bessell (1979) 
relations. Parenthetically, we note that a more careful scrutiny of the photometric data than in 
CAL4 resulted in a better definition of the period - color relation. Except for the short-period 
system RW Com, all systems are located within the band marked in Figure ||. 

The expected spread in the predicted values of My for CAL5, estimated from the Monte 
Carlo experiment described above, is shown in Figure |(| The continuous and broken lines give 
the spread in absolute magnitudes encompassing 68.3% and 95.4% of all cases; these intervals 
can be identified with one- and two-sigma ranges. As described above, the spread stems from 
the combined uncertainties in the coefficients and from the observed spread in (B — V). It is 
surprisingly small indicating that, in principle, the standard mean error of predictions of My 
should be at the level of 0.1. The observed mean error is larger, a = 0.22 (vertical bar in 
Figure |6|), indicating that the Monte Carlo estimated errors may give overly optimistic estimates. 
Of particular importance are three deviating systems with good parallaxes, V759 Cen, SW Lac 
and TY Men, which are clearly too bright than predicted by the calibration and which mainly 
contribute to the increase of the mean error a. We have checked the input data for these three 
systems and see no obvious reason why they deviate and what is causing the high luminosities, 
since they also do not give any clear hints of a third, hidden parameter dependence, as will be 
discussed in Section |||. 

4. Do we see a third parameter in the luminosity calibration? 

The new luminosity - period - color relation may be accurate enough to reveal systematic 
deviations of calibrating systems with properties deviating from the general trend. It was stated 
in CAL1 that the mass-ratio q would be the most likely third parameter in the calibration. Thus, 
we may expect that the calibration should be fully valid for systems with an "average value" of 
q, and that systems with mass-ratios deviating noticeably from the average value would show 
systematic residuals AMy from the calibrating plane. In addition to q, the fill-out factor / and 
the orbital inclination i may influence the luminosity to a certain extent. The first two, q and /, 
are related to structural properties and are known (Maceroni and van't Veer 1996) to correlate 
with the spectral type and total mass (early type systems have usually small mass ratios and large 
fill-out factors), so that we might expect that their contributions are already absorbed in the color 
term. 

We have collected the system parameters (q, f, i) for the majority of the systems used in 
CAL5. The data are given in Table [l], details on the source of data are given in the Appendix. The 
residuals AMy from the calibration are plotted versus each of these parameters in Figure 0. The 
formal linear fits, giving coefficients as in Table ||, indicate no significant dependencies for q and /; 
this table gives also results of bootstrap re-sampling experiments. The lack of systematic trends 
is also obvious in Figure [?]. Only the i dependence shows a 0.15 mag variation in the interval 60° 
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to 90° which is well covered by observed systems, indicating the possibility of a gravity-induced 
brightening at the poles. The slope per degree, +0.0053lo!oo4i> ^ s significant at the level slightly 
higher than one-sigma. Since this effect is of the order of the accuracy of the calibration, it can be 
neglected in statistical studies. 

In addition to the geometrical parameters (q, f, i) one could consider metallicity as a hidden 
unknown of the calibration. Unfortunately, nothing is known about differences in metallicities 
within the sample of the local W UMa-type systems. Similarly, usage of any correlated indicators, 
such as spatial velocities, would not be advisable as one would expect a statistical rather than a 
deterministic relation between deviations in My and the values of the space velocity vectors. 



5. The poor thermal contact and semi-detached systems 

As was described in CAL1 and CAL4, the existence of a period - color relation causes an 
insufficient separation of the color and period dependencies so that the results crucially depend 
on the data for "outlying" systems which are contact systems with unusually long periods for 
their colors. Such systems appear close to the red-color (broken line) boundary in Figure |2[ 
They are important as they can stabilize the solutions by lifting the period - color degeneracy. 
For the previous calibrations in CAL1 to CAL4, the system V371 Cep (NGC 188 V5) was 
used in this special role. The current Hipparcos sample defined by eMy < 0.5 includes a wide 
range of combinations of periods and colors, so that it is no longer necessary to rely on extreme 
systems such as V371 Cep. In fact, we are for the first time in a position to check if the 
calibration is indeed valid for such systems. As indicated by the OGLE data (Rucinski 1997b), the 
poor-thermal-contact (PTC) and possibly semi-detached (sd) systems mimicking good contact are 
very rare in space, constituting only some 2% of all contact systems. They are very important for 
our understanding of pre- and in-contact evolution. We used V371 Cep in the previous versions of 
the luminosity calibration to improve stability of solutions, but it has never been proven that the 
luminosity calibrations of the CAL series are applicable to such systems at all! . 

There are very few binaries with reliable values of My among the PTC/sd systems. In 
addition to V371 Cep, one can consider here RT Scl which seems to be a genuine semi-detached 
system (Hilditch and King 1986, Banks et al. 1990), perhaps also the red system RW Com 
(Milone et al. 1987) which is included in our calibration, although with very low weight. If 
we apply the new calibration to such systems, their absolute visual magnitudes (derived from 
spectroscopic and photometric analysis) are usually fainter than predicted. The deviations amount 
to 0.80 magnitude for V371 Cep, 0.40 magnitude for RT Scl, while RW Com deviates by as much 
as 0.95 magnitude. 

The sample of the PTC/sd systems is clearly too small to draw definite conclusions, but 
we point out that low luminosities (for a given color) are actually expected for such systems. 
By comparing mixing of colors of Main Sequence stars (clearly the extreme case of no thermal 
contact) with colors of genuine contact systems, one sees that the PTC/sd systems should be 
always bluer than systems in good geometrical and thermal contact (Figure |8|). The observed 
colors will obviously depend on how good or poor the thermal contact actually is, so that their 
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values cannot be predicted easily. But the result is clear: The luminosities estimated from the 
calibration are expected to be always too high, which seems to be confirmed by results described 
above. It would be of great importance to confirm this prediction on the basis of a larger sample 
of PTC/sd systems. 

6. Conclusions 

We have presented a luminosity calibration of contact binaries based on Hipparcos parallaxes, 
which covers a wide range in periods and colors. It is based on the {B — V) color index and 
supersedes previous calibrations based on that color index. It nominally predicts absolute 
magnitudes to an accuracy of ±0.1, although deviations of some systems are larger, indicating a 
wide, non-Gaussian distribution of uncertainties. For the current data, the weighted mean error is 
a = 0.22. The calibration can be written as: 

My = -4.44 log P + 3.02(5 - V) Q + 0.12 (5-a) 

We have investigated the influence of mass-ratio, fill-out factor and inclination, and have found 
that the first two parameters are, within the accuracy of the calibration, completely absorbed 
by the color term. Only the inclination, which is independent of the physical configuration of a 
system, shows a noticeable effect, as it is expected from the gravity brightening. Nevertheless, this 
effect is so small that it can be neglected in statistical studies. Close binary systems appearing 
in the period domain of the W UMa-type systems, but showing unequally deep eclipses (either 
poor-thermal contact or semi-detached systems) appear too faint relative to the calibration by 
roughly 0.5 magnitude, but the sample is much too small to draw definite conclusions. 

The accuracy of the present calibration is no longer limited by the parallax data but, 
paradoxically, by the lack of reliable photometric data. It should be noted that an uncertainty in 
the color index (B — V) of 0.03 contributes a deviation of 0.1 magnitude in the predicted My. 

HWD gratefully acknowledges support from the sabbatical visitors program at Space 
Telescope Science Institute and SMR acknowledges a research grant from the Natural Sciences and 
Engineering Council of Canada. The research has made use of the SIMBAD database, operated 
at the CDS, Strasbourg, France. 

Appendix: Sources of input data 

In the following, information is given on the source of apparent magnitudes at maximum, 
color indices and system parameters of the W UMa systems used in CAL5. 

AB And. V magnitude and (B — V) index from Landolt (1969). He quotes (B — V) = 0.90 as 
the color index at total eclipse. Using the color index of the comparison, and differential colors 
of the contact system, we get 0.865 for the (fairly poorly determined) maximum and 0.90 for the 
minimum. The (b — y)-index suggests {B — V) = 0.88, which was adopted for CAL5. The system 
parameters were taken from the detailed study of Hrivnak (1988), 
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S Ant. V and {B — V) from Popper (1957), system parameters from Russo et al. (1982). 

V535 Ara. V and (B — V) from Schoffel (1979), system parameters are averages of the values 
given by Schoffel (1979) and Leung and Schneider (1978). 

TZ Boo. This system shows a very unstable light curve. V and (B — V) averaged from data 
by Hoffmann (1978), mass ratio given by McLean and Hilditch (1983), inclination estimated by 
Al-Naimiy and Jabbar (1987). 

AC Boo. (B — V) given by Mauder (1964). Several B, V light curves are available, but none 
refers to a comparison star with known B, V magnitudes. The V magnitude at maximum was 
calculated from the Hp magnitude. Mass ratio from Hrivnak (1993), inclination and fill-out factor 
from Mancuso et al. (1978). 

CK Boo. Asian and Derman (1986) measured UBV light curves relative to HD 128128. 
Stromgren photometry of HD 128128 (Olsen 1994) yields V = 7.897, b-y = 0.348, from which 
(B — V) = 0.55 was derived. Mass ratio from Hrivnak (1993). 

RR Cen. V and (B — V) from Chambliss (1971), system parameters averaged from Mochnacki 
and Doughty (1972) and King and Hilditch (1984). 

V752 Cen. V and (B - V) from Sistero and Castore de Sistero (1973), also (B - V) from 
Duerbeck (1997). Averaged color index used here. System parameters taken from Barone et al. 
(1993). 

V757 Cen. V and (B - V) from Cerruti and Sistero (1982), also (B - V) from Duerbeck (1997). 
Averaged color index used here. System parameters taken from Maceroni, Milano and Russo 
(1984). 

V759 Cen. (B — V) was derived from Bond's (1970) b — y value. Maximum V magnitude was 
calculated from the Hp magnitude. 

VW Cep. Linnell (1982) gives B, V light curves relative to BD +74° 889, whose B, V magnitudes 
were taken from Sorgsepp and Albo (1974). System parameters were taken from Hill (1989) and 
Mochnacki (1981). 

RS Col. V and (B — V) from McFarlane and Hilditch (1987), who also determined system 
parameters. 

RW Com. V and (B — V) from Milone et al. (1980), system parameters from Milone, Wilson 
and Hrivnak (1987). 

e CrA. V and {B — V) from Tapia (1969), mass ratio from Goecking and Duerbeck (1993), 
inclination and fill-out factor by Tapia and Whelan (1975). 

SX Crv. V and (B — V) are averaged from Sanwal et al. (1974) and Scaltriti and Busso (1984). 

V1073 Cyg. The (B - V) index is given by Mendoza, Gomez, & Gonzalez (1978). Several light 
curves are available, but none refers to a comparison star with known B, V magnitudes. The V 
magnitude at maximum was calculated from the Hp magnitude. The system parameters were 
taken from the study of Ahn, Hill and Khalesseh (1992). 
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RW Dor. V and (B — V) from Marton, Grieco and Sistero (1989), who also give system 
parameters. The parameters used here are average values of the previous study and that of 
Kaluzny and Caillault (1989). 

AP Dor. V magnitude at maximum estimated from Eggen's (1980) observations, (B — V) taken 
from Przybylski and Kennedy (1965). 

BV Dra and BW Dra. V and (B — V) averaged from Yamasaki (1979) and Rucinski and 
Kaluzny (1982). The system parameters were taken from Kaluzny and Rucinski (1986). 

YY Eri. V at maximum calculated from the Hp magnitude. (B — V) taken from Eggen (1967). 
System parameters were derived by Nesci et al. (1986). 

SW Lac. V and {B — V) from Rucinski (1968). System parameters were derived by Leung, Zhai 
and Zhang (1984). 

XY Leo. V from Hrivnak (1985). The {B — V) index is an average of the determinations of Koch 
(1974), Hrivnak (1985), and Duerbeck (1997). System parameters were derived by Hrivnak (1985) 
and Barden (1987). The latter found that XY Leo is a quadruple system. The V magnitude has 
been increased by 0.05 to account for the M-type pair. 

AP Leo. V at maximum calculated from the Hp magnitude. (B — V) taken from Duerbeck 
(1997). System parameters were derived by Zhang, Zhang and Zhai (1992). 

UV Lyn. V and (B — V) from Bossen (1973). A more recent derivation of system parameters 
was carried out by Markworth and Michaels (1982). 

TY Men. Naqvi and Gr0nbech (1976) obtained uvby light curves. V magnitude averaged from 
both maxima; (B — V) calculated from average (b — y) index. The system parameters were derived 
by Lapasset (1980). 

UZ Oct. V and (B — V) taken from Castore de Sistero, Sistero and Candellero (1979). The 
system parameters were derived by Lapasset and Sistero (1984). 

V502 Oph. V and (B — V) taken from Wilson (1967). System parameters are averaged from 
results of King and Hilditch (1984) and Maceroni et al. (1982). 

V566 Oph. Lafta and Grainger (1985) give B, V photometry relative to BD +4° 3558, whose 
UBVRI photometry is given by Castelaz et al. (1991). System parameters are averaged from 
results of Van Hamme and Wilson (1985) and Mochnacki and Doughty (1972). 

V839 Oph. Several B, V light curves are available, but none refers to a comparison star with 
known B, V magnitudes. V at maximum was derived from the Hp magnitude. (B — V) is 
an average of the value published by Koch (1974) and that calculated from the (b — y) index 
determined by Rucinski and Kaluzny (1981). The inclination is taken from an analysis by Niarchos 
(1989). 

MW Pav. V and (B — V) from Lapasset (1977), system parameters taken from Lapasset (1980). 

U Peg. No light curve relative to a comparison star with known B, V magnitudes exists. V at 
maximum was derived from the Hp magnitude. (B — V) was taken from Eggen (1967), system 
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parameters from the study of Zhai and Lu (1988). 

AE Phe. V and (B — V) from Walter and Duerbeck (1988), system parameters were averaged 
from the results of Van Hamme and Wilson (1985) and Niarchos and Duerbeck (1991). 

VZ Psc. V at maximum was extrapolated from Bradstreet's (1985) fragmentary UBVRI 
observations. Several measurements of (B — V) are available; a "median" value is that measured 
by Ryan (1989). System parameters were determined by Hrivnak, Guinan and Lu (1995). 

AQ Psc. V and {B — V) from Sarma and Radhakrishnan (1982). 

HI Pup. V at maximum was calculated from the Hp magnitude, (B — V) was taken from 
Duerbeck (1997). 

V781 Tau. V and (B — V) from Cereda et al. (1988). System parameters were determined by 
Lu (1993). 

W UMa. We used the B, V light curves by Breinhorst (1971), which were measured relative to 
BD +56° 1399, for which Oja (1984) gives UBV magnitudes. System parameters were taken from 
the study of Linnell (1991). 

AW UMa. V and (B — V) from Paczynski (1964). System parameters were averaged from the 
studies by Rucinski (1992) and Hrivnak (1982). 

GR Vir. V and {B - V) from Cereda et al. (1988). 
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(B - V) 

Fig. 1. — The color - magnitude relation for the Hipparcos systems with absolute-magnitude 
errors smaller than 0.25 magnitude (filled circles) and within 0.25 to 0.5 magnitude (open circles). 
The same symbols are used in all subsequent figures. The line gives the main sequence for single 
solar-type stars, with the Sun marked by a cross. The short vertical vectors point at locations 
predicted by the new calibration presented in this paper called CAL5. 
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Fig. 2. — The period - color diagram for the W UMa-type systems observed by Hipparcos (filled 
and open circles, as in Figure 1). The solid line is the blue-short-period envelope (BSPE) from 
CAL3, transformed to the (B — V) color. The broken line gives the red edge of the band which was 
used in a Monte Carlo experiment to estimate the spread in My ■ Its location was roughly estimated 
to include evolved systems with red colors (for a given period). The red edge was constructed by 
adding a linearly-increasing value to the BSPE, so that the band would widen from 0.3 at 0.2 days 
to 0.6 at 1.0 day. The BSPE is poorly defined for periods longer than about 0.6 - 0.7 day. 
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Fig. 3. — The period - color - luminosity relation for the Hipparcos systems. The inclined surface 
gives the new calibration CAL5. The mean weighted deviation of the observed points from the 
calibration plane is a = 0.22. The data projected onto the horizontal period - color plane yield the 
same period - color relation as in Figure 2. 




Fig. 4. — Deviations from the new calibration CAL5 are shown here as functions of the two 
principal parameters. The vertical error bars have lengths equal to twice the absolute-magnitude 
errors, as estimated from the uncertainties in the Hipparcos parallaxes. 
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Fig. 5. — Relations between the coefficients of the new calibration CAL5 established by a bootstrap 
re-sampling experiment. The naming of the coefficients is the same as in CAL1, in that BV in 
parentheses indicates that this is a (B — V)-based calibration. The solid contours encompass 68.3% 
of all bootstrap solutions of the coefficients, a level which is normally associated with the one- 
sigma standard-error uncertainty. The dotted lines give similar levels for 50%, which is normally 
associated with the probable errors. Note that separate "islands" in mutual correlations between 
the coefficients - which were present in the previous calibrations - no longer exist indicating higher 
stability of the present solutions. 
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Fig. 6. — The spread in the predicted values of My obtained from a Monte Carlo experiment for 
CAL5 is shown here by curves. It was derived by combining the spread in the calibration coefficients, 
as shown in Figure 5, with an assumed spread in colors. In practice, 10,000 calibration coefficients 
obtained in the bootstrap experiments were applied to 100 uniformly distributed random values 
of the colors (within ranges shown in Figure 2) for each value of the period in steps of 0.01 day. 
The solid lines give the one-sigma range (68.3% of all cases), whereas the broken lines give the 
two-sigma range (95.4% of all cases). The median for the random samplings is shown by the dotted 
line. The two vertical broken lines delineate the region defined by systems with the extreme values 
of the orbital period in the Hipparcos sample. The vertical bar gives the weighted mean deviation 
in the sample, a = 0.22. The filled and open circles as well as the vertical bars have the same 
meanings as in other figures. 
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Fig. 7. — Deviations from the calibration are shown here versus the physical parameters of the 
contact system: the mass ratio, q = M2/M1 < 1, the fill-out parameter, /, defined through 
equipotentials with / = for the inner common equipotential, and the orbital inclination, i, in 
degrees. The filled and open circles as well as the vertical bars have the same meanings as in 
previous figures. The dashed lines indicate the weighted least squares fits to the data. Note that 
(q, /, i) data are not available for all systems studied here. 
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Fig. 8. — Expected changes in V and (B — V) (upper panel) and V and (V — I) (lower panel) for 
good-thermal contact systems with the varying mass-ratio, q, are shown by the continuous lines 
with marks every 0.1 in q. The broken lines show simple mixing of colors for Main Sequence stars 
which can be taken as an extreme case of poor thermal contact. The calculations presented here 
are for the color of the primary component (B — V) = 0.75 and (V — I) = 0.8 and assume the 
slope of the Main Sequence My oc 4.85(5 — V) and My oc 4.35(V — I). The color changes for good 
thermal-contact systems have been computed following the analysis of Mochnacki (1981) assuming 
the nuclear mass - luminosity relation L oc M +4 ' 4 . Note that the largest color changes are observed 
for 0.5 < q < 0.6, and that they are accompanied by very small changes in brightness. Apparently, 
for such mass-ratios, the contact systems have secondaries contributing large radiating areas, but 
almost no luminosity. For q —> the systems resemble single stars, whereas for q — > 1 the stars again 
influence each other to a lesser degree as they become identical, with the combined Ay — > —0.75. 
The poor thermal contact systems are expected to be located between the continuous and broken 
curves in the figure. 
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Table 1. Contact binary stars observed by Hipparcos with eMy < 0.5 
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1 


.0 


0.40 


1.65 


+0.05 


RW Dor 


0.2855 


0.89 


0.00 


10.90 


11.07 


0.68 


0.12 


77.0 


8.9 


2 


.0 


0.49 


5.65 


+0.43 


AP Dor 


0.4272 


0.44 


0.00 


9.26 


9.37 








4.7 





.9 


0.42 


2.62 


-0.46 


BV Dra 


0.3501 


0.56 


0.01 


7.89 


7.94 


0.41 


0.11 


76.3 


14.9 


2 


.6 


0.38 


3.72 


-0.08 


BW Dra 


0.2922 


0.63 


0.01 


8.61 




0.28 


0.14 


74.4 


14.9 


2 


.6 


0.38 


4.44 


+0.08 


YY Eri 


0.3215 


0.66 


0.00 


8.16 


8.30 


0.40 


0.15 


82.5 


18.0 


1 


.2 


0.14 


4.44 


+0.14 


SW Lac 


0.3207 


0.75 


0.01 


8.54 


8.80 


0.88 


0.21 


79.0 


12.3 


1 


.3 


0.23 


3.96 


-0.58 


XY Leo 


0.2841 


0.99 


0.02 


9.55 


9.64 


0.50 


0.06 


65.8 


15.9 


1 


.8 


0.25 


5.49 


+0.02 


AP Leo 


0.4304 


0.62 


0.00 


9.30 


9.43 


0.30 


0.24 


79.9 


8.3 


1 


.7 


0.44 


3.90 


+0.28 


UV Lyn 


0.4150 


0.67 


0.00 


9.42 


9.58 


0.50 




68.0 


8.2 


1 


.6 


0.42 


3.99 


+0.15 


TY Men 


0.4617 


0.29 


0.03 


8.11 


8.17 


0.22 


0.22 


79.5 


5.9 





.6 


0.22 


1.87 


-0.52 


UZ Oct 


1.1494 


0.54 


0.03 


9.03 


9.14 


0.28 


0.26 


78.5 


3.8 





.8 


0.46 


1.84 


+0.45 


V502 Oph 


0.4534 


0.64 


0.01 


8.34 


8.51 


0.37 


0.23 


71.3 


11.8 


1 


.2 


0.22 


3.67 


+0.13 


V566 Oph 


0.4096 


0.43 


0.01 


7.45 


7.55 


0.24 


0.25 


80.3 


14.0 


1 


.1 


0.17 


3.15 


+0.04 


V839 Oph 


0.4090 


0.63 


0.09 


8.82 


8.96 






75.0 


8.1 


1 


.4 


0.38 


3.08 


-0.39 


MW Pav 


0.7950 


0.36 


0.05 


8.53 


8.71 


0.18 


0.46 


85.0 


4.8 


1 


.1 


0.50 


1.78 


+0.29 


U Peg 


0.3748 


0.62 


0.02 


9.47 


9.60 


0.33 


0.10 


76.0 


7.2 


1 


.4 


0.42 


3.69 


-0.13 


AE Phe 


0.3624 


0.64 


0.00 


7.56 


7.69 


0.39 


0.20 


86.3 


20.5 





.8 


0.08 


4.12 


+0.11 


VZ Psc 


0.2612 


1.14 


0.00 


10.15 


10.33 


0.80 


0.15 


48.0 


16.8 


2 


.1 


0.27 


6.28 


+0.13 


AQ Psc 


0.4756 


0.50 


0.03 


8.60 


8.66 








8.0 


1 


.3 


0.35 


3.02 


+0.05 


HI Pup 


0.4326 


0.60 


0.07 


10.33 


10.46 








5.4 


1 


.2 


0.48 


3.77 


+0.44 


V781 Tau 


0.3449 


0.58 


0.05 


8.55 


8.68 


0.41 


0.27 


66.5 


12.3 


1 


.4 


0.25 


3.84 


+0.08 


W UMa 


0.3336 


0.62 


0.00 


7.76 


7.85 


0.49 


0.10 


86.0 


20.2 


1 


.1 


0.12 


4.29 


+0.18 


AW UMa 


0.4387 


0.35 


0.00 


6.84 


6.91 


0.07 


0.65 


79.1 


15.1 





.9 


0.13 


2.73 


-0.03 


GR Vir 


0.3470 


0.57 


0.00 


7.80 


7.93 








18.8 


1 


.2 


0.14 


4.17 


+0.29 
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Table 2. Calibration Coefficients of CAL5 





a P(W) 


&BV 


a 0(BV) 


Average 


-4.44 


+3.02 


+0.12 


Median 


-4.42 


+3.08 


+0.10 


— 1 a from median 


-0.81 


-0.70 


-0.23 


+ 1 a from median 


+0.75 


+0.51 


+0.13 



Table 3. Influence of a third parameter on the luminosity (in magnitudes per unit value) 



Parameter 


Q 

zero-point 


slope 


/ 

zero-point 


slope 


i 

zero-point 


slope 


Average 


-0.014 


+0.05 


+0.037 


-0.05 


-0.35 


+0.0048 


Median 


-0.013 


+0.04 


+0.040 


-0.06 


-0.39 


+0.0053 


— 1 a from median 


-0.062 


-0.21 


-0.051 


-0.11 


-0.19 


-0.0041 


+ 1 a from median 


+0.061 


+0.21 


+0.050 


+0.15 


+0.30 


+0.0027 



